Search for: All records

Gene duplication is a fundamental evolutionary mechanism that contributes to biological complexity and diversity [6]. Traditionally, research has focused on the duplication of gene sequences [23]. However, evidence suggests that the duplication of regulatory elements may also play a significant role in the evolution of genomic functions [8, 21]. In this work the evolution of regulatory relationships belonging to gene-specific-substructures in a GRN are modeled. In the model, a network grows from an initial configuration by repeatedly choosing a random gene to duplicate. The likelihood that the regulatory relationships associated with the selected gene are retained through duplication is determined by a vector of probabilities. That is to say that each gene family has its own probability of retaining regulatory relationships. Occurrences of gene-family-specific substructures are counted under the gene duplication model. In this work gene-family-specific substructures are referred to as subnetwork motifs. These subnetwork motifs are motivated by network motifs which are patterns of interconnections that recur more often in a specialized network than in a random network [15]. Subnetwork motifs differ from network motifs in the way that subnetwork motifs are instances of gene-family-specific substructures while network motifs are isomorphic substructures. These subnetwork motifs are counted under Full and Partial Duplication, which differ in the way in which regulation relationships are inherited. Full duplication occurs when all regulatory links are inherited at each duplication step, and Partial Duplication occurs when regulation inheritance varies at each duplication step. Note that Full Duplication is just a special case of Partial Duplication. Moments for the number of occurrences of subnetwork motifs are determined in each model. In the end, the results presented offer a method for discovering gene-family-specific substructures that are significant in a GRN under gene duplication. 

We examine the collective behavior of single cells in microbial systems to provide insights into the origin of the biological clock. Microfluidics has opened a window onto how single cells can synchronize their behavior. Four hypotheses are proposed to explain the origin of the clock from the synchronized behavior of single cells. These hypotheses depend on the presence or absence of a communication mechanism between the clocks in single cells and the presence or absence of a stochastic component in the clock mechanism. To test these models, we integrate physical models for the behavior of the clocks in single cells or filaments with new approaches to measuring clocks in single cells. As an example, we provide evidence for a quorum-sensing signal both with microfluidics experiments on single cells and with continuousin vivometabolism NMR (CIVM-NMR). We also provide evidence for the stochastic component in clocks of single cells. Throughout this study, ensemble methods from statistical physics are used to characterize the clock at both the single-cell level and the macroscopic scale of 10⁶cells. 

Estimating transmission rates is a challenging yet essential aspect of comprehending and controlling the spread of infectious diseases. There are various methods available for this purpose, each with its own assumptions, data requirements, and limitations. This paper introduces a phylogenetic approach called transRate, designed to estimate inter-population transmission rates. The phylogenetic method, which maintains statistical consistency under the multi-population Susceptible-Infected-Recovered (SIR) model, integrates genetic information with traditional epidemiological approaches. This integration improves the accuracy of transmission rate estimates, facilitating more effective disease control and prevention strategies. Simulation analyses validate the precision of transRate in estimating transmission rates. With the growing abundance of public databases for genomic sequences, transRate is becoming more prevalent in tracking and preventing the spread of such diseases. 

Abstract Utilizing a microfluidic chip with serpentine channels, we inoculated the chip with an agar plug withNeurospora crassamycelium and successfully captured individual hyphae in channels. For the first time, we report the presence of an autonomous clock in hyphae. Fluorescence of a mCherry reporter gene driven by aclock-controlled gene-2 promoter(ccg-2p) was measured simultaneously along hyphae every half an hour for at least 6 days. We entrained single hyphae to light over a wide range of day lengths, including 6,12, 24, and 36 h days. Hyphae tracked in individual serpentine channels were highly synchronized (K = 0.60-0.78). Furthermore, hyphae also displayed temperature compensation properties, where the oscillation period was stable over a physiological range of temperatures from 24 °C to 30 °C (Q₁₀ = 1.00-1.10). A Clock Tube Model developed could mimic hyphal growth observed in the serpentine chip and provides a mechanism for the stable banding patterns seen in race tubes at the macroscopic scale and synchronization through molecules riding the growth wave in the device. 

Abstract The presence of Arbuscular Mycorrhizal Fungi (AMF) in vascular land plant roots is one of the most ancient of symbioses supporting nitrogen and phosphorus exchange for photosynthetically derived carbon. Here we provide a multi-scale modeling approach to predict AMF colonization of a worldwide crop from a Recombinant Inbred Line (RIL) population derived fromSorghum bicolorandS. propinquum. The high-throughput phenotyping methods of fungal structures here rely on a Mask Region-based Convolutional Neural Network (Mask R-CNN) in computer vision for pixel-wise fungal structure segmentations and mixed linear models to explore the relations of AMF colonization, root niche, and fungal structure allocation. Models proposed capture over 95% of the variation in AMF colonization as a function of root niche and relative abundance of fungal structures in each plant. Arbuscule allocation is a significant predictor of AMF colonization among sibling plants. Arbuscules and extraradical hyphae implicated in nutrient exchange predict highest AMF colonization in the top root section. Our work demonstrates that deep learning can be used by the community for the high-throughput phenotyping of AMF in plant roots. Mixed linear modeling provides a framework for testing hypotheses about AMF colonization phenotypes as a function of root niche and fungal structure allocations. 

System biology relies on holistic biomolecule measurements, and untangling biochemical networks requires time-series metabolomics profiling. With current metabolomic approaches, time-series measurements can be taken for hundreds of metabolic features, which decode underlying metabolic regulation. Such a metabolomic dataset is untargeted with most features unannotated and inaccessible to statistical analysis and computational modeling. The high dimensionality of the metabolic space also causes mechanistic modeling to be rather cumbersome computationally. We implemented a faster exploratory workflow to visualize and extract chemical and biochemical dependencies. Time-series metabolic features (about 300 for each dataset) were extracted by Ridge Tracking-based Extract (RTExtract) on measurements from continuous in vivo monitoring of metabolism by NMR (CIVM-NMR) in Neurospora crassa under different conditions. The metabolic profiles were then smoothed and projected into lower dimensions, enabling a comparison of metabolic trends in the cultures. Next, we expanded incomplete metabolite annotation using a correlation network. Lastly, we uncovered meaningful metabolic clusters by estimating dependencies between smoothed metabolic profiles. We thus sidestepped the processes of time-consuming mechanistic modeling, difficult global optimization, and labor-intensive annotation. Multiple clusters guided insights into central energy metabolism and membrane synthesis. Dense connections with glucose 1-phosphate indicated its central position in metabolism in N . crassa . Our approach was benchmarked on simulated random network dynamics and provides a novel exploratory approach to analyzing high-dimensional metabolic dynamics. 

Abstract A consensus species tree is reconstructed from 11 gene trees for human, bat, and pangolin beta coronaviruses from samples taken early in the pandemic (prior to April 1, 2020). Using coalescent theory, the shallow (short branches relative to the hosts) consensus species tree provides evidence of recent gene flow events between bat and pangolin beta coronaviruses predating the zoonotic transfer to humans. The consensus species tree was also used to reconstruct the ancestral sequence of human SARS-CoV-2, which was 2 nucleotides different from the Wuhan sequence. The time to most recent common ancestor was estimated to be Dec 8, 2019 with a bat origin. Some human, bat, and pangolin coronavirus lineages found in China are phylogenetically distinct, a rare example of a class II phylogeography pattern (Avise et al. in Ann Rev Eco Syst 18:489–422, 1987). The consensus species tree is a product of evolutionary factors, providing evidence of repeated zoonotic transfers between bat and pangolin as a reservoir for future zoonotic transfers to humans. 

Abstract We determined the macroscopic limit for phase synchronization of cellular clocks in an artificial tissue created by a “big chamber” microfluidic device to be about 150,000 cells or less. The dimensions of the microfluidic chamber allowed us to calculate an upper limit on the radius of a hypothesized quorum sensing signal molecule of 13.05 nm using a diffusion approximation for signal travel within the device. The use of a second microwell microfluidic device allowed the refinement of the macroscopic limit to a cell density of 2166 cells per fixed area of the device for phase synchronization. The measurement of averages over single cell trajectories in the microwell device supported a deterministic quorum sensing model identified by ensemble methods for clock phase synchronization. A strong inference framework was used to test the communication mechanism in phase synchronization of quorum sensing versus cell-to-cell contact, suggesting support for quorum sensing. Further evidence came from showing phase synchronization was density-dependent. 

The use of a microwell microfluidic device allows separating single cells and tracking single cells data. The measurement of single cell fluorescent intensity trajectories in the microwell device supported a deterministic quorum sensing model identified by ensemble methods for clock phase synchronization. A strong inference framework was used to test the communication mechanism in phase synchronization of quorum sensing versus cell-to-cell contact, and the results lent support for quorum sensing.